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VARIATIONAL ENERGY PRINCIPLE FOR 
COMPRESSIBLE, BAROCLINIC FLOW 
I. FIRST AND SECOND VARIATIONS 
OF TOTAL KINETIC ACTION 
Lawrence A. Schmid 

Goddard Space Flight Center, Greenbclt, Maryland 20771 

ABSTRACT 

As a preliminary to the construction of the effective total free-energy functional for 
general inviscid flow that is carried out in the following paper, the special case of a cold 
gas in the absence of external force fields is considered. Since the only energy involved 
is kinetic energy, the total kinetic action (i.e., the space-time integral of the kinetic energy 
density) should serve as the total free-energy functional in this case, and as such should be 
a local minimum for all possible fluctuations about stable flow. This conjecture is tested 
by calculating explicit, manifestly covariant expressions for the first and second variations 
of the total kinetic action in the context of Lagrangian kinematics. The first variation 
vanishes if the motion of every particle is rectilinear, and the second variation depends on 
the velocity shear. If the Lagrange surfaces are parameterized so as to correspond to in- 
compressible flow, a generalization of Kelvin’s Energy Theorem results, in the sense that 
the total kinetic action is shown to be a minimum not only for potential flow, but also 
for steady Beltrami flow and for unsteady flow that satisfies Euler’s equation with spatially 
homogeneous Bernoulli constant. Which of these cases arises depends on the imposed 
boundary and terminal conditions, which are more general than those of Kelvin’s Theorem. 
Finally, the general question of the correlation between physical stability and the convex- 
ity of any action integral that can be interpreted as the total free-energy functional of the 
flow is discussed and illustrated for the cases of rectilinear and rotating shearing flows. 
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The type of instability in question is characterized by an evolution that is consistent with 
the preservation of a specified average-velocity profile for the flow, and so corresponds to 
the onset of ordered turbulent motion within a given large-scale flow. 
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VARIATIONAL ENERGY PRINCIPLE FOR 
COMPRESSIBLE, BAROCLINIC FLOW 
I. FIRST AND SECOND VARIATIONS 
OF TOTAL KINETIC ACTION 

L INTRODUCTION AND SUMMARY 

It is the purpose of this paper and the one that follows (which will be referred to as 
Paper II) to construct a variational principle for inviscid, compressible, baroclinic flow in 
a gravitational field that can serve as a basis for direct solution of the kinds of problems 
encountered in dynamic meteorology. By “direct solution” is meant the trial-and-error 
procedure in which the desired solution is expressed in terms of parameterized trial func- 
tions, and the numerical values of the parameters are determined by minimizing the action 
integral of the variational principle. 

In Paper II this objective will be accomplished by constructing an action integral 
whose first variation vanishes when the known equations of motion are satisfied, and whose 
second variation has the form of the time-average change in the total free-energy of the 
system that is produced by arbitrary fluctuations in the different variables of the problem 
about the extremal * trial functions that satisfy the equations of motion. Thus S# can be 
interpreted as the time-average total free-energy of the flow, and its construction amounts 
to the extension of the formalism of thermodynamic potential functions to include a fluid 
in arbitrary motion. A stable flow corresponds to a local minimum (in parameter space) 
of this total free-energy, whereas an unstable flow corresponds to a saddle point whose 
properties are Correlated to the physical nature of the instability. 

The various thermodynamical questions will be left to Paper II. The present paper 
deals with the kinematics of the problem, and the only energy discussed will be the kinetic 
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energy. The kinetic action integral ^ is the time-average of the total kinetic energy, and 
the variational principle based on W can be regarded as the fluid generalization of the 
Least Action Principle for particle motion. It is shown in Section IV A that the flow for 
which ^ is an extremum is pure rectilinear motion of the individual particles. This cor- 
responds to the motion of a cold compressible gas in the absence of any external field. 

In Section V it is shown that, if the trial functions are parameterized so that the 
fluctuations are density-preserving, then the extremal flow is no longer limited to recti- 
linear motion. Rather, it is the flow that satisfies Euler’s equation/for spatially uniform 
Bernoulli constant, a condition that can result from irreversible turbulent mixing. In other 
words, the extremal flow is of the kind that results after mixing has caused the fluid to 
“forget” its initial energy distribution. For time-independence, the extremum corresponds 
to Beltrami flow, which includes potential flow as a special case. For density-preserving 

fluctuations, the extremum is always a minimum. Thus the variational principle based on 
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'W' with density-preserving fluctuations is a generalization of Kelvin’s Energy Theorem, ’ 
which says that the total kinetic energy is less for potential flow than for any other in- 
compressible flow having the same normal velocity component at the boundary. The 
reason that the extremal flow discussed in Section V is more general than the steady po- 
tential flow of Kelvin’s Theorem is that the imposed boundary conditions are correspond- 
ingly more general. 

Section VI A is a discussion of the relation between physical stability and the second 
variation of any action integral which, like the one discussed in Paper II, can be in- 
terpreted as the total free-energy of the flow. (The action integral ‘W, whose extremal 
flows are restricted by the complete absence of any potential energy, is a rather uninter- 
esting special case.) It is shown that the requirement that >0 for all possible 
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fluctuations about the extremal flow constitutes a sufficient, but not necessary, condition 
for stability. The explicit expression for 5 (2) W that is given in Section IV A is applied 
in Section VI B to derive sufficient criteria for stability of parallel shearing and rotating 
flows. 

It is pointed out in Sections IV C, D and VI A that any stability criterion derived 
from a variational energy principle refers to a different kind of instability from that which 
is predicted by a normal-mode analysis of the linearized inviscid differential equations of 
motion. This latter variety, which is called “laminar instability” in this paper, is typified 
either by the onset of a self-excited growing sinuous oscillation, or else by the transition 
to an entirely different mode, such as the onset of convection. In both cases, the original 
velocity profile of the base flow is destroyed by the instability. In the case of the criterion 
S* 2 * Stf > 0, however, the instability, which is called “internal instability,” is constrained 
to preserve the average-velocity profile of the base flow. The jet streams in the upper 
troposphere provide an example of both types of instability. A large-scale undulation, 
which completely changes the original velocity field of the jet is an example of “laminar 
instability,” whereas the high-level clear-air turbulence that is encountered near the jet 
maximum, and consists of eddies embedded in the jet wind and carried with it, is an ex- 
ample of “internal instability.” In Section VI B 2, it is noted that observations of high- 
level clear-air turbulence confirm that, for this type of instability, the internal stability 
criterion is indeed the appropriate one to use. In Section IV D it is pointed out that both 
kinds of instability are necessary for a complete description of turbulent flow. 

The fact that a stability analysis based on a variational energy principle involves dif- 
ferent modes of instability from those predicted by an analysis of the differential equations 
of motion is intimately related to the fact that, in order to solve the differential equations, 
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it is necessary to specify the initial velocity of every particle, whereas in the case of a vari- 
ational principle, it is effectively the time-average velocity that must be specified. This 
means, as explained in Section IV B, that a variational analysis corresponds to an eigen- 
value problem which in general does not have a single unique solution as in the case of 
the differential-equation analysis, but rather has a whole ensemble of dynamically accept- 
able solutions. This difference, and its physical interpretation, is elaborated in Sections 
IV C and D. 

At the end of Section VI A, it is pointed out that one way of partially bridging the 
gap between the extreme idealizations represented by the “laminar” and “internal” sta- 
bility criteria is to relax the terminal and boundary conditions imposed on a variational 
analysis. A simple example involving the vertical mixing of different layers in a horizontal 
wind shear is given in Section VI C. 


In Section VII it is shown that the first variation of W is invariant under a Galilean 
transformation, but not the second variation. This means that the convexity of W (or of 
any action integral of which it is a part) can be interpreted as a stability criterion only if 
the fluid velocity is referred to an appropriately chosen Galilean frame which must be 
identified with the “momentum reservoir” with which the fluid system under study is in 
equilibrium. The interpretation of the terminal and boundary conditions of the variational 
principle in terms of a momentum reservoir in interaction with the flow system under study 
is made in Section IV D. 


The central analytical result of this paper is the explicit expression for 5 (1 + 


- < 


that is given in Equation (48) of Section IV A. As a preliminary to deriving this expres- 
sion, the formalism of Lagrangian kinematics is developed in Section II. It is absolutely 
necessary that the Lagrange, rather than the Euler, description be used. The physical 
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reason for this is that, in order to find a functional that behaves like the total free-energy 
of the flow, it is necessary that the fluctuations of the trial functions imitate the physical 
fluctuations that actually occur in an inviscid fluid. First and foremost, this means that 
the fluctuations must be such that the conservation equation is identically satisfied for all 
trial functions, not just for the set of extremal functions. This is possible in the Lagrange 
formalism, which uses the three families of Lagrange surfaces A A (<r, t), (A =1, 2, 3), as 
the basic kinematical variables, but it is not possible in the Euler description, which treats 
the four functions n(&, t) and Vl(#, t), (j = 1, 2, 3), as primitive and independent. In the 
Lagrange formalism n and V are secondary functions that are expressed in terms of the 
derivatives of A A in such a way that the conservation equation is identically satisfied for 
arbitrary A A (#, t). 

In Section III the necessary functional differentials are derived. It is emphasized that, 
in order to obtain second-order functional differentials that can be assigned an absolute 
physical meaning in the sense of a stability criterion, the expressions must be manifestly 
covariant with respect to transformations both of the inertial coordinate system and of the 
system of Lagrange surfaces. The derivation can be done more straight-forwardly in an 
inertial frame than in the convected frame. However, in applying the variational prin- 
ciple to the numerical solution of a problem, only an evaluation of the total action integral 
is involved (rather than the derivation of a local functional differential), and this can be 
done as conveniently in the convected frame as in an inertial frame (much more so, in 
fact, if a free boundary is involved). 


II. KINEMATICS OF LAGRANGE DESCRIPTION 
A. Reciprocal Tripod Fields 

The entire kinematics of the Lagrange description is generated either by the three 
functions A a (#, t), or by the three functions XJ(A a , t) (or the single vector function 
X(A, t». The respective spatial derivatives of these functions generate two reciprocal tri- 
pod fields G a and G A ; 

G A = ei(3A A /axj) t = V A a ; (1) 

G a = OX/3A A ) t S 8 a X . (2) 

These tripod fields satisfy the following identities: 

G a • G b = 5»; G a G a = C*G A = 7 (3a, b) 

where the summation convention with respect to the index A is used, and S A and 7 are 
respectively the Kronecker delta and the idem dyadic, The tripod field G A can be ex- 
pressed directly in terms of the vectors G B , and hence in terms of derivatives of A a ( X , t), 
rather than those of X(A, t), by means of the following relation: 

G a = G b X G c /j£ (4) 

where is the functional determinant (Jacobian) whose elements are the derivatives 
G A s 3 A a jdx J. A similar relation expresses the G A tripod directly in terms of the de- 
rivatives 9X3/3 A A . 


These tripod fields suffice to generate all of the mathematical quantities encountered 
in curvilinear tensor analysis.^ For example, the metric tensor g AB and the Christoffel 
symbol T^ B are given by 


g ab = g a 


G 


B »■ 


f ab = (3a G b ) • GC = O a 9 b X) • G c . 


(5) 

( 6 ) 
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Thus the option exists either to work directly with the tripod fields, or else to use G A B 
and T AC in the conventional Riemannian formalism. In general, the first course is the 
simpler. 


If U is an arbitrary vector, its contravariant form U A in the convected frame is con- 
structed by projection onto the G A tripod, and the covariant form U A by projection onto 
the G a tripod : 

U A = u • G A ; U A = U • G a • (7a, b) 

Although any vector can always be cast into either contravariant or covariant forms, it 
turns out that every physical vector has a preferred or “natural” form in terms of which 
the equations of motion assume their simplest form. For example, in the case of fluid 
dynamics the natural form of the velocity vector is V A whereas the natural form of the 
molar momentum vector is MV A . It is typical that canonically conjugate pairs such as 
these have natural forms that arel sontragredient so that a direct contraction yields a scalar. 

IS 

The definition of the Lagrangian surfaces A a at the initial time t r is arbitrary and so 
also is the resulting Jacobian (J^) r Thus the normalization of the Lagrangian surfaces 
must be specified. In this paper and in Paper II they will be normalized so that at tj 
the volume contained in every parallelepiped defined by neighboring pairs of A A -surfaces 

A 

that differ by unity is equal to the molar volume V of the fluid. Because these surfaces 

move with the fluid, this normalization will be maintained for all later time. Because 
X . 

V" 1 = n, the mole density, it follows that for the specified normalization 


■ /s x / 9XA ■ 

v = I * Hdet W’ (8a) 

n = }} = (VAl) • (VA2 x V.U) (8b) 

(where the indices have been underlined in order to indicate that they are numerical values 
of A rather than of j). 
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The following identities are very useful: 


dV/9G A = VG a ; 9n/9G A = nG A ; 

(9a, b) 

9 A (VG A ) = 0; V- (nG A ) = 0; 

(10a, b) 

9G a /9G b = -G b G a ; 9 G a /9G b = -G B G A ; 

(Ha, b) 

D t G A = G a • V V; D t G A = -(VV) • G A 

(12a, b) 


where D t is the material time derivative (defined in (13) and (45c) below). 

B. Convected Frame 

In the convected frame of the fluid the independent coordinates are (A a , t) and the 
kinematics is generated by the vector function X(A a , t) which gives the time-dependent 
position vector X (referred to an inertial frame) of the particle that is identified by the 
Lagrange parameters A A . The molar volume V is given in terms of X(A, t) by (8a), and 
the velocity V is given by 

V = O t X) A = D t X. (13) 

Because of the normalization chosen for the A-surfaces, the volume element d 3 (A) in the 
convected frame can be identified with the infinitesimal mole number dN. 


Thus the total instantaneous kinetic energy W is given by 


W 


= / V2M(b t X) A • (B t X) A d 3 (A) = f WdN (14) 


'(A) •'N 

where M is the molar mass, and W = VM\ • V is the molar kinetic energy, and N is the 
total number of moles of fluid contained in the region of integration. The total kinetic 
action W is defined by 


w = f F wdt = f 


y 2 MO t X) A V O t X) A dt d 3 (A) 


(15) 


(At, N) 


where At = t F - tj is the time range of integration. Note that except for the factor At, 
can be interpreted as the time average of the total kinetic energy: 



W/At, 


( 16 ) 


In fact, the time unit can be chosen so that At = 1, in which case 9(f and < W> t are numer- 
ically equal, and the two can be regarded as conceptual equivalents. 


C. Inertial Frame 

In the inertial (i.e., laboratory) frame the independent coordinates are (X, t). (Note 
that x is used to denote the position vector when it serves as the aggregate of independent 
spatial coordinates whereas X is used to denote the dependent variable X(A, t) that gives 
the position vector of a given particle. The two are, of course, numerically equal, i.e., 
x - X.) The position vector X will be taken to be expressed in Cartesian coordinates so 
that the volume element d 3 (x) is equal to the physical differential volume element d V. 


The kinematics is generated by the three families of Lagrange surfaces A A (x, t). The 
mole density n is given in terms of VA a by (8b). As indicated in (7), the velocity V can 
be written in the form 

V = V a G a ; V a - V • G A . (17) 

Because V, and hence V A , is not a primitive variable (as in the Euler description) it 
is necessary to express V A in terms of A A (X, t). Such a relation results from substitu- 
tion of dA A = 0 and dX = V dt into dA A = G A • dX + (3 t A A ) x dt. Making 
use of (17), the result is ,r . 

V A = -O t A A ) x ; V = -O t A A ) x G a . (18) 

The vectors G A are expressed in terms of VA® by means of (4). Thus the expression to 
be used fcr V in Lagrangian kinematics in an inertial frame is 
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3 

V = - £ O t A A ) x ( VA ® x VA c )/n; (A, B, C cyclic) 


(19) 


A = 1 


where n is given by (8b). When n and V are given by (8b) and (19) respectively, the con- 
servation equation is identically satisfied: 

(3 t n) x + V • (nV) = 0 . (identity) (20) 

This means that, for an arbitrary choice for the three functions A A (<r, t), the Lagrangian 
kinematics will always satisfy conservation of matter. This is not the case for the Eulerian 
description in which t) and n(dc, t) are primitive and independent variables. This is 
the reason why fluctuations in the functional forms of A a (#, t) may be regarded as gen- 
erating valid descriptions of the physical fluctuations in a turbulent fluid, whereas fluctu- 
ations in the functional forms of V(x, t) and n(iC, t) may not be so regarded. 


In the inertial frame, the total kinetic action W' is given by 
<9* 


-L 


VisMV • Vnd 3 (x)dt = / wd 3 (x)dt 
(v. At) •'(V, At) 


( 21 ) 


where V and n are given by (19) and (8b) respectively, and V is the total instantaneous 
volume of integration, and w = nW is the instantaneous kinetic energy density. 

D. Steady Flow 

A limited class of steady flows can be treated very easily in the Lagrangian formalism 
by eliminating the time-dependent functions A A (#, t) in favor of the time-independent 
functions X A 0r) by means of the relations 

A a (<*, t) = X A (ir) - V A t (22) 

where the V A are three given constants. Then 

= J A = = (VX 1 ) • (7X- x VX-), (23) 

G A = (VA A ) t = VX A = g A (24) 
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( 25 ) 


G. = G A x G B /n = g A x g B /n == g , (A, B. C cyclic) 

)) * 

V=-O t A*v4 = vA « A IJ 

3 

= £ V a (VA b X VA c )/n. (A, B, C cyclic) 

A = l 

It follows directly from (26) and the constancy of the V A that 

V • (nV) = 0 . (identity) 


(26) 


(27) 


From (23) and (26) it is evident that the total kinetic energy W is 

W = J VM\ • Vnd 3 (x) = / 1 / 2 M(g AB V A V B ) j£ d 3 (x) , 

w 


(28) 


where g AB = g A • g B and g A is given by (25). 


The relation of the A-formalism to the usual stream function formalism can be shown 
by considering the special case V A = (0, 0, 1). Then (26) becomes 

nV = VAl x V\2 (29) 

which is the standard way of expressing the flux density nV in terms of two stream func- 
tions A - (x) and A- (x). The A-formalism differs, however, from the stream function 

1 ! 2 • 3 

formalism in that the former is based on the three functions A- (x), A- Or), and A- (x) 
whereas the latter is based on the three functions A - (x), A - ( x ), and n(x) (where n is set 
equal to a constant for an incompressible fluid). The fact that the stream function for- 
malism suppresses one of the three X A (x) in favor of n(x) means that the complete map- 
ping necessary to transform from the a; -frame to the A-frame is not available, so the stream 
function formalism can only be applied in the -frame whereas the A-formalism can be 
applied in both the * -frame and the A-frame. 
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In the X-frame (28) becomes 

w= X ,/jMvAvB ©-E) d3(x) - 

For the special case V A = (0, 0, 1) this becomes 


(30) 


3 

and it is evident that X- has taken over the role played by t in the time-dependent 
formalism. 


(31) 


Only a limited class of steady flows can be represented in the form (22). It will now 
be shown that the physical nature of this limitation consists of the fact that for closed 
(i.e., re-entrant) steady flows represented by (22) the average of the covariant strain-rate 
tensor taken over any closed streamline is zero. The same is true of many open-ended 
flows. This means that only steady flows that are shearless in this average sense can be 
represented by (22). The proof starts with the observation that the strain-rate tensor can 
be referred to the “natural” coordinate system defined by the X-surfaces as follows: 

ej(9iv k + 9 k VJ)e k = g A (V A V B + V B V A )g fi (32) 

where V A V B is the covariant gradient of the velocity defined by 


vAvB = ^ (i^ + v °) 


(33) 


where r£ D is the Christoffel symbol defined on the tripods g A and g A in the manner 
analogous to that given in (6). For the flow under consideration 3V B /9X C = 0 and ap- 
plication of Gauss’ theorem leads to the following result: 


/ 


(V A V B +V B V A )nd 3 (x) 


■l 


g AB nV • d^. 


(34) 
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When the volume and surface integrals involved are taken over the volume and surface of 
a very fine elemental flow tube, the desired proof results. For simplicity, it will be stated 
for the special case V A = (0, 0, V-). Then (34) becomes: 

' I 

1 f X l 3 

<V A V B + V B V A > 3 = — — I (V A V B + V B V A )dX- 

x3 (X| - Xf ) J X 3 (35) 

= -v^(g£ B - g f B )/(x| - \-) 

where X~ and X- are respectively the values of X- at which the fluid enters and leaves 
the elemental flow tube. Thus the tube-average of the covariant strain-rate tensor is pro- 
portional to the difference between the values of the metric tensor g A B at the two ends 
of the tube. If the flow is closed, this difference vanishes, and so too does the tube-average 
of the strain-rate tensor. The same is true for open-ended flow for which the geometry 
defined by the X-surfaces is the same at the entry and exit surfaces of the flow. In this 
sense, we can say that (22) describes steady flow that on the average is shearless. 

When the flow is independent of the X- -coordinate, the covariant strain-rate tensor 
vanishes at every point in the flow. This is the case for axially-symmetric rotational flow 
in which X- is the azimuth angle. In this case V- is the angular velocity, and the fact that 
it is a constant means that the fluid moves as if in rigid rotation. Thus in this case (22) 
is unable to represent those non-rigid vortices, such as the potential vortex, that involve 
non-vanishing velocity shear. 

It is possible to generalize the representation given in (22) so that the V A are no 
longer constant and hence restricted to shearless flows. In fact, besides admitting arbitrary 
shear, both V A and X A can be made time-dependent. This formalism, which combines 

the, advantages of the Euler and Lagrange representations, will not be developed in this 

■ ■ ■.-■■■■' '■r-VW ■ r 

h 

paper. 
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It should be emphasized that the general time-dependent A-formalism, which will be 
employed throughout the remainder of this paper, does not share any of the restrictions 
discussed above in connection with (22). The general A-formalism is perfectly capable of 
describing either steady or time-dependent shearing flow. However, steady shearing flow 
in the A-formalism is characterized by the fact that one or more of the G A fields is time- 
dependent, even though the velocity field as seen in the inertial frame is time-independent. 


ra. FIRST AND SECOND FUNCTIONAL DIFFERENTIALS 


A. Calculational Methods 

If the objective were merely to calculate 5 (1 + 2) W , i.e., the first and second vari- 
ations of the total kinetic action, then this could be done either in the convected frame 
using the expression for W given in (15), or in the inertial frame using (21). In either 
case the result would be expressions that would show the sensitivity of the total kinetic 
action to changes in the functional dependence of the three functions that describe the 
kinematics. 

If, however,- the objective is to obtain a local second-order expression that can be 
interpreted as a stability criterion, then the desired second-order expression must have 
an absolute significance that is independent of the choice of coordinates and of the defi- 
nition of the A a . This means that both the independent coordinates and the displace- 
ments (i.e., variations) appearing in the expression must be referred to an inertial frame. 

The reason for this is that second-order virtual changes in kinetic energy are being used 
as a measure of stability, and kinetic energy must be defined with respect to an inertial 
frame. Correspondingly, the Prinicple of Virtual work can be related to changes in kinetic 
energy only if the displacements involved are the particle displacements 5X. Then the 
virtual work done on the particle can be related to its change in kinetic energy. Thus it 
is necessary to derive an expression involving 5X(#, t) whose coefficients are also functions 
of Or, t). 


In the convected frame the second of these requirements is violated because the in- 
dependent variables are (A a , t), and in the inertial frame the first requirement is violated 
because the variations involved are 5 A a (it, t) rather than 5X(;c, t). Thus in either case 
it would be necessary to transform the original expression for the variation in order to 
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give it absolute significance. It turns out to be more straight-forward to perform the cal- 
culation in the inertial frame, and then eliminate 6A A (ac, t) in favor of 8X(x, t). It can 
be shown that these two displacements are related by 

-5X(<r,t) = G a 8A a + ViG c 5 A a 5A b + 0 3 (5A) (36) 

where G A and T AB are calculated from A B Gr, t) by means of (4) and (6) respectively. 

The overhead bar is used to indicate the particular arbitrarily chosen functional dependence 
of A B (ir, t) about which the expansion is carried out, Thus 5 A a = A A («, t) - A a (#, t) 
where the functional form of A a O r, t) differs slightly from that of A a (jjc, t). The minus 
sign on the left side of (36) results from the fact that, if the fluid is carrying particles 
toward a region of larger A a - values, then the A a -values that flow past an observer fixed 
in the inertial frame appear to decrease in value. 

Except for the minus sign, the first-order part of (36) is just the statement of the 
standard transformation law of a vector. The extra second-order contribution results from 
the fact that -5X and 5 A a are not really different expressions for the same vector, but 
rather they are two physically different vectors. The vector SX is the displacement of a 
single particle as it traverses many Cartesian coordinate planes, whereas the 5 A a are the 
changes in Lagrange parameters of many different particles that flow past an observer fixed 
at the intersection of a single given set of Cartesian planes. 

It is evident from (36) that an expression that is of first order in SX will appear to 
be of both first and second order in 8 A a , and correspondingly for a first order expression 
in 5 A a that is converted into one involving SX. Thus, when the first-order variation of 
W expressed in terms of 5 A a is converted into an expression in SX, this expression will 
contain second-order terms which will contribute to the stability criterion. 
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Note that none of these subtleties enter into a first-order calculation of the variation, 
or a direct application of the Euler-Lagrange equations^ to find the differential equations 
that must be satisfied if the first variation of the action integral is to vanish. What little 
attention the second-variation has received in the classical mathematical literature has con- 
centrated on the problem (Jacobi’s Accessory Problem^) of establishing a sufficient con- 

. I 

dition for a positive-definite second variation (in the form of the existence of solutions 
to a certain set of partial differential equations that can be derived from the action inte- 
gral). From this point of view, the question of the absolute significance of the variation 
never arises. 

It should be noted that, except for the sign on the left side of (36), this relation is 
identical to the one used in Riemannian geometry to introduce geodesic coordinates. ^ In 
the present context, the introduction of geodesic coordinates amounts to redefining the 
Lagrange parameters A a in such a way that at a given instant t 0 (the moment at which 
the variations are being calculated) the redefined Lagrange surfaces coincide everywhere 
with the Cartesian planes of the inertial frame. This coincidence will not hold for t =£ t Q , 
but at t = t Q the Christoffel symbol for the geodesic coordinates (i.e., the redefined 
Lagrange parameters) is everywhere zero. This fact can be used to simplify the calculation 
of + 2) ^ a t t Q , and because t Q is arbitrarily chosen, the expressions that are de- 
rived can be applied at any instant of time. 



The results given below in Section III B were calculated by the geodesic method just 
outlined, and in addition by an entirely independent, somewhat less abstract method. This 
second method involves using (36) directly to replace 5 A a by 5X. It also involves taking 
into account that the tripod vectors G A and G A fluctuate when the functional forms of 
A A (<r, t) fluctuate. In constructing an absolute vector, the components of the vector 
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after the change in the functional forms of A a Or, t) must be contracted with the appro- 
priate tripod vectors as calculated also after the change. This change in tripod vectors is 
most easily expressed in terms of 5X by making use of several simple identities. If X and 
X = X + 5X are the old and new particle positions respectively, and (G A , G A ) and (G A , 


G A ) are the old and new tripods, then 

G a s (3X/3A a ) = 3(X - 5X)/3 A a = G a - 3 a 5X . (37) 

Contracting this with G A and using (36), the desired identity results: 

G a G a = 7 - 75X (38) 

which shows that G A G A can be regarded as a generalization of the idem dyadic 7 - 
G a G a . The reciprocal of (38) gives another identity: 

G a G a = (G b G b ) _1 = 7 + 75X + (V5X) • (V5X) + 0 3 (fiX). (39) 

Similar identities can be found that involve 5 A a instead of 5X: 

G A G A = 3T+(7SA a )G a (40a) 

* 7 + (75A a )G a + (76A a )(3 a 5A b )G b + 0 3 (5A) ; (40b) 

G A G A = 7 - (75A a )G a + (75A a )(3 a 5A b )G b + 0 3 (5A) (41a) 

= 7 - (75A a )G a (41b) 

where 3 A = G A • 7and 3 A = G A • 7. Expressions for 6 (1 + 2 ^G A and d*- 1 + 2) G A re- 
sult from the contraction of (39) and (38) with G A and G A respectively: 

■fid + 2 ) G a W G a - G a = G a • 76X + G a • (75X) • (75X) (42a) 

fid + 2 >G a .= G a - G a = -(75X) • G A ....... . (42b) 

As an example of the transformation of an expression from one involving 6 A a into 


one involving 5X, consider the expression for fid + 2 > S(A a , t) = S(A a , t) - S(A a , t) 
where S(A A , t) is a specified functional relation. In this case it is easy to show that it is 
necessary to use nothing more than (6) and (36) in order to change the expression for 
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6 (1 + 2) S from one involving 8 A a into one involving 8X: 

6< 1 + 2 > S(A, t) = 8 A x (d^S) + 1 / 2 5A A 8A B (6 A a B S) (43a) 

= “5X • VS + y 2 (5X5X : VVS). (43b) 

B. Functional Differentials of n, V, and w 

The functional differential of the mole density n results from substitution of A a = 
A a + 8 A a into n = det (8 A a /3x j ) and subsequent conversion of the resulting expression 
in 5 A a into one involving 5X. This was done by both of the methods outlined above, 
with the following results: 

5 (1 >n = _ V. (n 5X) ; (44a) 

5(1 + 2) n = _ V* 8 (1 + 2 >n (44b) 

where 

5 (1 + 2 >n = n(5X + 6X • V5X) - %V • (5Xn5X) (44c) 

or 

8 (1 + 2 >n = (ii + V&5< x >n) (8X + V 28 X • V8X) . (44d) 

From (44b) it follows that 5 (1 + 2 * n may be regarded as being caused by a density flux 
5U + 2) n produced by the field of particle displacements 8X. 

The functional differential of V results from substitution of A A = A a + 5A a into 
(19) and subsequent conversion of the resulting expression in 5 A a into one involving 5X. 
This was done by both the geodesic coordinate method, and by the more direct method 
involving the identities given in (37)-(42), and the following results were obtained: 

5d) v = D t 5X - 5X • VV ; (45a) 

5(1 + 2 ) y = D t (5X + 8X • V5X) - 8X • V(V + D t 5X) (45b) 

where . 

D t = (d t ) x + V • V. S 


Since V = D t X, it follows that D t 6X is just the change in V that is required by the change 
5X. The term -SX • VV in (45a) can be interpreted as the advected change in V as 
follows: The matter that is located at the observation point x after the displacement SX 
was originally located at the point x - SX before the displacement. If velocity behaved 
like an embedded property of the matter, then -5X • V V would be the change in veloc- 
ity that would result simply because the matter at the observation point had been changed. 

The expression for 5 (1 + 2) w results simply from substitution ofn = n + 5 (1 + 2 >n 
and V = V + 6 (1 + 2) V into w = VinMV • V where S* 1 + 2) n and 5 (1 + 2) V are given 
by (44b) and (45b) respectively: 

5(1 + 2) w - _ v . [5(1 + 2) n( ^ _ i/ 2 5 X . VW + MV • D t 8X)] 

+ nMV • D t (8X + 6X • V5X) + %nM(D t 8X) • (D t SX) (46) 

- y 2 nM(5X5X :VVV) • V. 

All of the terms on the right have simple intuitive interpretations. The first term says that 
the effective flux of kinetic energy density is just equal to the density flux 5 (1 + 2) n 
multiplied by the molar kinetic energy W, corrected for the average advected molar kinetic 
energy -V 28 X • VW and the increase MV • D t 5X that results from the velocity increase 
D t 5X. The second and third terms on the right give the expected increase in kinetic 
energy that results from the velocity change caused by 5X. The last term on the right 
side of (46) can be interpreted in terms of advection as follows: If V were an embedded 
fluid property like the function S in (43), then the second-order change in it produced 
by 5X would be 5^ 2) V = V2(5X8X : VW) and the corresponding change in molar kinetic 
energy would be 5^ 2 -* W = HM(5X5X : V V V) • V. Thus the last three terms in (46) can 
be interpreted as the mole density n multiplied by the increase in molar kinetic energy at 
the observation point that exceeds the second-order energy that is advected to the ob- 
servation point. 
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IV. VARIATIONAL ENERGY PRINCIPLE BASED ON TOTAL KINETIC ACTION 
A. General Expression for + °H' 

The expression for 5*- 1 + 2 ^ ^ follows directly from (21) and (46). By means of 
partial integration the expression can be cast into the form presented below. This will be 
written in terms of [S^ 1 + [S^ 1 + 2)< ^] s , and [5 (1 + 2 ^^¥'] v which are respec- 

tively the terminal, surface, and volume contributions, the last of these being represented 
as the sum of five contributions; 


4 

[5< 1 + 2 Wj v = [6 (1 >^] v + £ [S (2 >3¥1^ N) . (47) 

N = 1 

The general expression for 8 (1 + 2) 3¥', which in no way restricts the choice of trial func- 
tions (indicated by an overhead bar) about which the fluctuations take place, is the 
following: 


-L 


5(1 + 2)^ = I [6 (1 + 2) W ] d 3( x) dt 

f (v,At) 


(48a) 


*T 

= [5< 1 + 2 >^] t + [5(i +2 W] s + [ 6 ^ W] y + X! [ 5<2)< ^v N) 


N = 1 


where 


[5 (1 


i 2 >W) T = / [nMV*(5X + 8X« V8X)l t = tp d 3 (x) 
j In MV • (5X + 5X • V5X)] t = tj d 3 (x) ; 


(48b) 


[ 6 (i + 2 )<^] s = _ / (W - V45X • VW + MV • D t 8X + MA • 5X) 5 (1 + 2 >n • d^dt ; 

(&* At > (48c) 
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[5 (1 W] V = 


[5< 2 )^](i) = 


[6 


[ 5 (2)^](3) _ 


[S< 2 ^]^ 4) = 


i 


■'(v , At) 


I 


(v,At) 


L 


(v , At) 


l 


(v , At) 


-X 


(v , At) 


MA • 5Xnd 3 (x)dt; 


MK« (1) n)A • 5X + V2n(5X8X : VA)] d 3 (x)dt ; 


y 2 M(D t 5X) • (D t 5X) n d 3 (x) dt ; 


VSM(8X8X : VA) n d 3 (x) dt ; 


%M(8X8X VVV) • Vnd 3 fx)dt 


(48d) 


(48e) 


(48 f) 


(48g) 


(48h) 


where A — D t V - 3 t V + V • VV is the local acceleration. 

Because the displacement field SX(x, t) is arbitrary, it follows that a necessary con- 
dition for 5 < 1 ) Df = 0 is that 

r\j rv nu r^j 

A = D t V = (3 t V) x + V • VV = 0 (49) 

throughout the space-time volume (V, At). The tilde, which has replaced the overhead 
bar, indicates quantities that are appropriate to extremal flow, i.e., for which the condition 
6 (1 W = 0 is satisfied. The condition (49) is necessary for = 0, but is not suf- 

ficient because the first-order variation of the terminal and surface integrals given in (48b) 
and (48c) must also be made to vanish. These integrals will be discussed in Section IV B. 

The necessary condition for the first-order vanishing of 5 Dt' that is given in (49) is 
the Euler-Lagrange equation^" for a variational principle for which the Lagrangian density 
is just the kinetic energy density expressed in terms of Lagrangian kinematics. The extremal 
flow described by (49) is rectilinear in the sense that every particle trajectory in space-time 
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is a straight-line. Note, however, that this does not require that the trajectories be 
parallel. 


If (15) is written in the form 

= 



w dt 


dN 


( 50 ) 


it is evident that the variational principle based on W is just the fluid generalization of 
the Least Action Principle of Maupertuis and Euler. ^ The particle action is simply replaced 
by the molar action, and the total action a )Y is the sum of the individual molar actions. 


Such a variational principle can be given a geometric interpretation in the sense that 
it generates a geometry whose geodesics are determined by the Euler-Lagrange equation 
given in (49), which are just straight lines. The molar action /w dt plays the role of an 
arc-length. This interpretation is very well known** in the case of particle dynamics. This 
interpretation can be extended to a more general variational principle, such as the one de- 
veloped in Paper II, in which the kinetic action W" is replaced by a more general total 
action S/ . The Euler-Lagrange equations generated by S/ determine the system of “geo- 
descis” to be associated with S/. For a complete specification, however, it is necessary to 
specify the end-points Xj and X F of each particle trajectory (geodesic) as well as the 
Euler-Lagrange equation that it obeys. The specification of these end-points for all par- 
ticle trajectories that traverse the space-time region over which W or S/ is integrated 
constitutes the terminal conditions, and part of the boundary conditions, that must be 
imposed on the variational principle. Before discussing these in Section IV B, some ob- 
servations will be made regarding ] y . 

For the extremal trial functions that satisfy (49), it is obvious that the integrals given 
in (48e) and (48g) vanish, with the result that 
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[5 (2) ^] v = / y2M[(D t 5X) 2 - (5X5X:VVV)-V]nd 3 (x)dt, (51) 

mV, At) 

where the overhead tilde indicates quantities appropriate to the extremal trial functions 
for which (49) is satisfied. The first term in the integrand of (51) is positive-definite re- 
gardless of the properties of the extremal flow, but the sign of the second term depends 

ru 

on the spatial dependence of \(x, t), and it will be this term that determines the stability 
of parallel shearing flow, which will be discussed in Section VI Bl. The positive-definite 
expression (D t SX) z can be made small by specifying 5X(#, t) to have the weakest pos- 
' sible t-dependence. Because the end-points Xj and X F must be held fixed during the 
variation, a certain time-dependence of 5X(#, t) is unavoidable if the trajectory X(A, t) 
is displaced from its extremal value X(A, t). However, this minimal time-dependence can 
be made arbitrarily small by making the time interval At = t F - tj arbitrarily large. For 
this reason, in calculating the stability of shearing flow, the first term in (51) will be ne- 
glected. It will be seen that this is a conservative procedure in finding a sufficient con- 
dition for stable flow. 

Note that n d 3 (x) = d 3 (x) - d 3 (A) = dN is just the infinitesimal mole number. 
Thus (51) could be expressed as an integral over mole number (and time). 


B. Terminal and Boundary Conditions 

It was noted above that the Euler-Lagrange equations of any action integral could 
be regarded as the determining equations of a system of “geodesics” whose complete de- 
termination further requires that the end-points Xj and X F of every geodesic (i.e., par- 
ticle trajectory) be specified. If or is being discussed in the context of the con- 
vected frame this means that Xj(A) = X(A, tj) and X F (A) = X(A, t F ) are specified func- 
tions of A a , and all trial functions X(A, t) must be equal to these specified functions for 
t = t t or t = t F . In the context of an analysis in an inertial frame A^ (x) = A a (*, t t ) and 


Ap(#) = A A (ir, t F ) are specified functions of x, and all trial functions A a (jc, t) must 
be equal to these at tj and t F . In other words, 

[SA A l t = t[; , F 5 [A a (*, t) - A**, t)J, = , F = 0 . (52a) 

The 5X that appears in (48b) is determined by 5 A a by means of (36). Thus it follows 
from (52a) and (36) that 

SXj * 5 X f = 0 . (52b) 

It follows from (52) and (48b) that 

[ 5 U + 2)^] x = o (53) 

if all trial functions A a (#, t) reduce to the specified functions A a (x) and A a {x) for 
t = tj and t = tp respectively. Note that (53) is valid even if the V that appears in (48b) 
is not the extremal velocity field that satisfies the Euler-Lagrange equations. 

In discussing the surface integral given in (48c), it is necessary to distinguish between two 
kinds of surfaces - access surfaces (designated by ^), and bounding surfaces (designated by 
«^ B ). Access surfaces are those through which fluid passes. These are further subdivided into 
entry and exit surfaces depending on whether the fluid enters or leaves the integration volume 
V. Bounding surfaces, either rigid or free, are those through which fluid never passes. 

The surface condition to be imposed on an access surface is really just an extension 
of the terminal conditions. This is illustrated in Figure 1. Figure 1 A shows a case in 
which the space-time integration volume (V, At) has been chosen in such a way that no 
fluid enters or leaves V. Thus only bounding surfaces are involved. The end-points 
of all the trajectories he in the terminal volumes Vj and Vp. If, however, the space-time 
volume of integration were chosen as shown in Figure lb, then for the same set of tra- 
jectories many of the end-points would fall on rather than on either of the terminal 
volumes. It is still necessary that the end-points of the trajectories be specified, but now 
this requirement has been partially shifted from the terminal volumes to the access surface 


PARTICLE 

TRAJECTORIES 


PARTICLE 

TRAJECTORIES 


Figure 1. Space-Time Integration Volume. (A) Trajectories Cross Bounding Hypersurface Only at Terminal Volumes V, 
and V F but not at Bounding Surface ^g. (B) Trajectories Also Cross Bounding Hypersurface at Access Surface ^ A . 





. From the perspective of four-dimensional space-time, the terminal surfaces V x and 
V F and the totality of bounding and access surfaces constitute a single hypersurface that 
encloses the four-volume (V, At). The combined terminal and boundary conditions then 
require that the point at which any particle trajectory crosses the bounding hypersurface 
must be specified, and must be the same for all trial functions. Wherever a trajectory is 
tangent to the hypersurface, no condition is imposed except for the one that defines 
tangency: 

5X • = 5X • d^ B = 0; (5X x d^ B ) is arbitrary. (54a) 

Wherever the condition 5X • d«^ = 0 is not fulfilled, d«^ is by definition an element of 
an access surface and the requirement is 

SX = 0 if 5X • =£ 0 . (54b) 

The requirements stated in (54) together with (52b) constitute the conditions to be im- 
posed on the trial functions on the bounding hypersurface. 

It follows from (44d) that 

§U + 2 > n = 0 if SX = 0 . (55a) 

Moreover, it can be shown that 

5 (1 + 2) n • && = 0 if 5X • = 0 . (55b) 

From (48c), (54), and (55) it follows that 

[fid + 2 >4 ¥':} s = 0 (56) 

if the surface conditions (54) are satisfied. From (48a), (53), and (56) it is evident that 

4 

gd + = [5d)<^] v + J] [5(2)^](N> ( 57 ) 

N = 1 

if the combined hypersurface conditions of (52) and (54) are satisfied. This is true for 
fluctuations about any arbitrary set of trial functions. For fluctuations about the extremal 
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set, the first-order term on the right-hand side of (57) vanishes (as do two of the four c 
second-order terms). 

The surface conditions (55) guarantee that the total matter N contained within V at 
every time instant is the same for all trial functions, i.e. 

fid + 2 >n = 0 for t,<t<t F (58) 

if the surface conditions stated in (55) are satisfied. 

The terminal conditions (52) can be given an equally simple interpretation if it is 
noted that the time-average velocity <V> t of a particle is given by 

<V>t = (Xp - Xj) / (tp - tj). (59) 

(For trajectories traversing an access surface instead of a terminal volume, the <V) t j s given 
by an obvious modification of (59).) The point to be made is that the imposed hyper- 
surface conditions guarantee that, regardless of the form that the trial functions might 
assume, the time-average velocity of every particle is always the same. Thus 

5<V> t = 0 (60) 

for every given set of A A values. 

The fact that N and (V) t are held constant constitutes a form of “normalization” of 
the problem that has the effect of bounding W from below. Thus the “Fluid Least Action 
Principle” based on W with the hypersurface conditions discussed above, has many fea- 
tures in common with the more familiar variational principle that can be constructed from 
any linear eigenvalue equation.^ If H is an energy operator and E is one of its eigenvalues, 
and <p the corresponding eigenfunction, then the eigenvalue equation is H0 - E (p , and 
this can be solved by finding the stationary points of the functional <0l H l0>/<0 10> where 
the brackets represent appropriate integrals. The normalization <0l0> = 1 can be imposed 
as a constraint, in which case the action integral is simply <0 ! IT 1 0 >. This integral is the 
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analog of <9 ? , and the normalization condition on 0 is the analog of the hypersurface con- 
ditions (52) and (54). Both types of normalization bound the respective action integrals 
from below. Each is also necessary in order to remove an ambiguity: The normalization 
of 0 is not determined by the linear eigenvalue equation, and so must be specified; and the 
Euler-Lagrange equations of W' do not suffice for a unique specification of the trajectories 
until the end-points are specified. 

This analogy indicates that the variational principle based on <9?, or a generalization 
of it based on the total action may be regarded as a complicated (nonlinear) eigen- 
value problem. This way of thinking will now be used to make the point that, whereas 
the extremal flow for 'W with given hypersurface conditions is unique, this is not in 
general the case for a more complicated variational principle. 

C. Non-Uniqueness of Extremal Flow 

In the case of the variational principle based on <W, the fact that the extremal tra- 
jectories must be straight lines in space-time means that, once the end-points of these 
lines on the bounding hypersurface are specified, the solution is uniquely determined 
because there exist no other straight lines that can be passed through the same end-points. 

In the case of a generalized variational principle of the kind that will be developed in 
Paper II, the extremal trajectories are no longer straight lines, and so it can no longer be 
asserted that a specification of the end-points on the bounding hypersurface, together 
with the Euler-Lagrange equations, suffices for a unique determination of the flow. In 
fact, the analogy with the eigenvalue problem discussed above indicates that this is not 
the case because an eigenvalue problem in general possesses a whole spectrum of solutions, 
each one having its own characteristic symmetry properties. 
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The point can be made in terms more appropriate to fluid dynamics as follows: If, 
in addition to the kinetic energy, various forms of potential energy are also included in 

the dynamics described by the action integral then wave motion is possible. Thus, for 

: f ■ 

given end-points on the hypersurface, in addition to the “smoothest pcissible flow” that 
connects these end-points and also obeys the Euler-Lagrange equations, it is also possible 
to obtain a dynamically acceptable flow by superimposing on the smooth flow one or 
more modes of wave motion. Besides the multiplicity that results from such superposition, 
there is the nonlinear “flip-flop” type of multiplicity that can involve a change in the 
topology of the flow. An example of this kind of multiplicity is the sudden appearance 
or disappearance of a vortex street or a bubble of back-flow embedded in the main flow. 

From the point of view of using a variational principle to carry out a direct numerical 
solution of a flow problem, this multiplicity of solutions would be interpreted to mean 
that the variational problem had been incompletely specified. In addition to conditions 
imposed at the bounding hypersurface, it would be necessary to specify sufficient “internal 
boundary conditions” in order to make the extremum unique. This would mean specifying 
such things as the amplitude and mode of any wave motion to be superimposed on the 
smooth flow. 

From a physical point of view, however, it must be recognized that such supplementary 
internal boundary conditions are not in general supplied by nature, and the multiplicity of 
solutions must be regarded as constituting an ensemble of dynamically acceptable flows. 

The role of this ensemble in the overall description of turbulent flow will be discussed in 
Subsection D below. 

The multiplicity of solutions that in general characterizes a variational principle 
stands in striking contrast to the uniqueness of the solution of the initial-value problem 
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that is based on the same set of Euler-Lagrange equations. In the initial-value problem 
for a fluid, instead of specifying X r and X F , or what amounts to the same tiling, Xj and 
<V> t , we specify Xj and the initial velocity Vj. That is, the average velocity in the vari- 
ational principle is replaced by the initial velocity in the differential-equation approach in 
which the Euler-Lagrange equations are solved directly. This replacement suffices to 
make the solution of the differential equations unique. In other words, although the same 
amount of information is specified in the two cases, the fact that it involves average velocity 
in one case as compared with initial velocity in the other has remarkable consequences. 

This difference is illustrated in Figure 2. 

There is a simple physical explanation of this difference. In the initial-value problem, 
by specifying initial velocity and particle position, we are also specifying energy and mo- 
mentum distribution, as well as all the spatial symmetries and topological features that 
characterize the flow at tj. None of these things are specified, however, in the average- 
value (i.e., variational) approach. For example, modes of flow having very different total 
energies can satisfy the same conditions on X T and <V> t . Thus, in the average-value 
approach, we are effectively specifying an ensemble of flows that is characterized by given 
average conditions, whereas in the initial-value approach we are discussing the deterministic 
time-evolution of a completely specified single flow. 

D. Physical Interpretation of Terminal and Boundary Conditions 

The requirement that the average value of the velocity of every particle must be held 
constant can be explained in terms of a thermodynamic analogy. If a system is immersed 
in a heat bath, or, more exactly, an entropy reservoir, then the free exchange of entropy 
between system and reservoir maintains the temperature of the system (or, more exactly, 
the time-average of the temperature, since fluctuations are always present) equal to the 
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constant temperature of the reservoir. Temperature and entropy are conjugate variables, 
and the behavior of an entropy reservoir is typical of any thermodynamic reservoir: The 
extensive variable is freely exchanged between system and reservoir, and this keeps the 
average value of the corresponding intensive variable of the system equal to that of the 
reservoir. In the case of a moving system in contact with a momentum reservoir, this 
contact would keep the average value of the velocity of the system equal to the given 
velocity of the momentum reservoir. The same argument applies even if the reservoir 
velocity is a function of the spatial coordinates such as <V> t . Thus the fluid described by 
the variational principle may be regarded as a system in contact with a momentum reser- 
voir characterized by the spatially-dependent velocity <V> t . 

Within this system various dynamically acceptable modes of motion are possible, and 
these constitute an ensemble of “internal” fine-scale motions that are consistent with the 
imposed “mean” motion <V) t (which could be a slow function of time as well as a function 
of the spatial coordinates). This decomposition of the velocity field into a mean velocity 
and a fine-scale “internal” velocity corresponds exactly to the usual way of discussing 
turbulent flow in terms of a mean motion that satisfies the Euler equation (in the invisdd 
limit) to which a Reynolds force term has been added, and a superimposed fine-scale 
motion that determines the Reynolds force. This force is given by the divergence of the 
Reynolds stress tensor whose components are given by (minus) the local time averages 
<(AVj)(AV k )) t where AV = V - <V) t is the deviation of the instantaneous velocity from the 
time-average velocity. Obviously, this stress tensor is different for the different modes of 
internal motion. The most stable mode will dominate. If this is “smooth,” i.e., if the 
velocity deviation AV and its derivatives are small, then the Reynolds force is small and 
the mean flow satisfies an equation that is nearly identical to the Euler equation, and so 
the flow is called “laminar.” If, however, the smoothest possible internal mode is not the 


most stable one and the mode with the least free energy is considerably less smooth, then 
the Reynolds force is significant, and the equation of motion for the mean flow will differ 
significantly from the Euler equation. In this case the flow is called “turbulent.” In such 
a case the flow problem could, in principle, be solved by the following iterative procedure: 
(1) The mean flow is solved using Euler’s equation without the Reynolds force term. (2) 
Using the calculated mean flow to specify <V> t , an appropriate variational principle is used 
to find the most stable fine-scale flow (which may or may not be smooth) that is con- 
sistent with the given <V) t . (3) The Reynolds force field for this fine-scale flow is calcu- 
lated. (4) This force field is added to the equation of motion for the mean flow, and 
this is solved again in order to determine (V) t more accurately. Step (2) is then repeated, 
etc. 

This argument shows that whether the mean flow is laminar or turbulent is intimately 
related to whether or not the smoothest internal mode is stable. Sufficient conditions for 
the stability of several simple smooth flows will be given in Section VI. 

The above argument also shows that a complete solution of turbulent flow involves 
both the initial-value approach (differential equations) and the average-value approach 
(variational principle). Instabilities can enter into either the mean-flow calculation or the 
fine-scale calculation. The first type is the familiar kind of stability problem that is much 
discussed in the literature, usually in the context of a normal-mode analysis. The fine- 
scale or internal stability, however, involves an ensemble and is necessarily a statistical con- 
cept. Ideally, it should be treated in a statistical context that assigns an appropriate sta- 
tistical weight to every member of the ensemble. However, it can also be treated, albeit 
more crudely, in a thermodynamic context that uses the total free energy of each flow 
mode to determine the most stable mode of the ensemble which is then regarded as the 
unique solution of the internal flow problem. 


36 


In those problems in which the smoothest internal mode is stable, the Reynolds force 
is small and the first step of the, iteration procedure outlined above suffices by itself to 
give a good solution of the large-scale flow. There exists another type of problem for 
which the second step is sufficient by itself. Dynamic meteorology provides many ex- 
amples of this type in which the flow configuration of interest, such as a convection cell, 
is embedded in a much larger general flow whose velocity profile is given. 

The thermodynamic approach to the determination of the most stable mode of an 
ensemble presupposes a knowledge of the functional that plays the role of the total free 
energy of the flow. In the case of incompressible flow, it is evident that this role must be 
played by the total time-average kinetic energy The reason for this is that the in- 
compressibility constraint decouples the compressive energy of the fluid from the kinetic 
energy, so only the distribution of kinetic energy is available to determine the flow. This 
special case will be discussed in the next section, and it will be shown that W does indeed 
possess the main characteristic of a free-energy functional, namely that it is a minimum 
for flow that is observed in nature, i.e., flow that satisfies the dynamical equations of 
motion. 
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V. GENERALIZATION OF KELVIN S ENERGY THEOREM 

Kelvin’s Energy Theorem^’ states that, for incompressible flow of a fluid having a 
spatially uniform density, the total kinetic energy is less for potential flow than for any 
other flow having the same normal velocity field at the boundary. Kelvin actually proved 
this theorem for time-dependent flow that is driven by the motion of an impermeable 
boundary, but all the integrals involved in his proof are time-independent spatial integrals, 
so the proof is the same as one for steady flow. Thus, in terms of the present discussion, 
Kelvin’s proof effectively involves a (non-Lagrangian) steady-flow treatment of the action 
integral in which the terminal volumes Vj and V F have been replaced by access surfaces, 
and the specification of A A (x, t) on these surfaces is constrained to be consistent with 
potential flow. The consequence of this constraint is the fact that only potential flow 
makes the first variation of the action integral vanish. If the specification of A a (x, t) on 
the access surfaces were allowed to be completely arbitrary, then the class of extremal 
flows would be larger. It will be shown below that it becomes the class of steady Beltrami 
flows. H The treatment given below is still more general, in that the full time-dependence 
of the action integral is taken into account, with the result that the class of extremal flows 
is still further enlarged. The treatment is further generalized in that density stratification 
is admitted, although the fluctuations are still constrained to satisfy the time-independent 
conservation equation. 

The nature of the implied constraint on the boundary conditions that is involved in 
Kelvin’s proof can be clarified by constructing an extended ensemble of steady flows 
corresponding to all possible specifications of normal velocity on the access surfaces that 
are consistent with a given fixed total flux through the surfaces. Then it can be shown 
(proof not given below) that, of all the extremal (i.e., Beltrami) flows in the extended 
ensemble, potential flow has the least total kinetic energy. 

39 


^CEDING 


page blank 


Not FTTMm* 


The generalized Kelvir, Theorem is obtained by specializing the general expression for 
5(1 + 2 )<^ given in (48) to the case for which the trial functions are parameterized in such 
a way that the condition 

0 5d) n = 5< 2 >n = 0 (61) 

is identically satisfied. Without making any use of (61 ), the expression for 5 (1 +2) w given 
in (46) can be cast into the following form: 

6<1 + 2>W = - V- {5 a+2) n[2(W + y 2 5 x W) + MV* 5 t V+MA * 5X]} 

+ iiD t [MV • (5X + 5X • V6X)] - W5 (1+2) n 

- M(n + 5^)n) 5X • [O t V) x - VX(VXV)] (62a) 

+ Mn {5X5X : V [O t V) x - VX(VXV)]} 

+ ViMn [(5 t V) • (5 t V) + (5 X V) • (5 X V)] 

where 

6 t V = D t 5X; 6 X V s -5X-VV; 5 X W = - 5X-VW. (62b) 

The desired expression is obtained by integrating (62a) over (V, At), making use of Gauss’ 
Theorem to convert the first term on the right side of (62a) into a surface integral, and 
converting the second term into terminal integrals by integrating over t in the A-frame. 
Upon imposing the incompressibility condition (61), and the terminal and boundary con- 
ditions discussed in Section IV B, the following expressions result: 


■ ■[* t* > «Vl Sn _ 0 - 

- / M[(9 t V) x - VX(VXV)] • 6Xnd 3 (x)dt, 

(63a) 

•'(v.At) 


16 (2) W'1 s „„ 0 = 

/ ; M {5X5X : V[(d t V) x - VX(VXV)]} nd 3 (x)dt 

MV, At) 1 

r 

(63b) 

+ 

I HMKStV) • (5 t V) + (5 X V) • (5 X V)] nd 3 (x)dt. 



Mv,At> 



At a given point in space the condition (61) represents a constraint on the spatial deri- 
vatives of 5X, but not on 5X itself. Thus, a necessary and sufficient condition that 
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fid ^ = Ois 


(a t V) x -VX(VXV) = 0 (64) 

throughout the space-time volume (V,At). The overhead bar has been replaced by a tilde 
in order to designate the extremal flow. Substitution of (64) into (63) yields 

[6 (1) #] 5n=0 = 0; [5< 2 >#] 5n=0 >0. (65a, b) 

The fact that the second variation is positive-definite means that ^ is a minimum for the 
choice of trial functions A a (x, t) that causes (64) to be satisfied. It is obvious that steady 
potential flow satisfies (64). However, even for steady flow, (64) admits the larger class 
of Beltrami flows. ^ (Non-potential Beltrami flow is characterized by streamlines that 
twist like the strands of a rope.) 

The physical significance of time-dependent flows that satisfy (64) becomes apparent 
when it is noted that the Euler equation for a compressible (or incompressible) inviscid 
fluid moving in the presence of a gravitational field 0( as, t) can be written in the form 

M[(3 t V)^ - VXfVXV)] = - V(H + M0+W) + TVS (66) 

where T is the absolute temperature and H, W, and S are respectively the molar enthalpy, 
molar kinetic energy, and molar entropy. (In the special case in which the given mole density 
i\(x, t) has the constant value n 0 , the right side of (66) reduces to - VB 0 where B 0 = f 5 /n 0 + 
U<p + W is the Bernoulli constant for unstratified incompressible flow.) 

If S and the molar Bernoulli constant, B = H + M0 + W, are both spatially uniform, 
the right side of (66) vanishes, and (66) reduces to (64). Thus the most general flow for 
which is a minimum for density-preserving fluctuations of the trial functions is flow 
for which the characteristic thermodynamic quantities S and B have become spatially 
uniform. This can happen through turbulent mixing. Because this is an irreversible pro- 
cess, and because complete spatial homogeneity is this final stage of this process, the flow 
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that satisfies (64) must be regarded as more stable, in the sense of the Second Law of 
Thermodynamics, than any other flow that satisfies the same boundary conditions but for 
which B and S are not spatially uniform. Thus maximum stability in the sense of the 
Second Law is equivalent to the minimum value of W that is consistent with the given 
boundary and terminal conditions. This equivalence is made self-evident if W'/At is re- 
garded as the time-average of the total free energy of the system, i.e., the energy that is 
convertible into useable forms. In classical quasi-static thermodynamics, the statement 
that stable equilibrium is characterized by a minimum value of the free energy of the sys- 
tem is one way of stating the Second Law. When time-dependence is admitted, it is 
evident that the instantaneous quantities of quasi-static thermodynamics must be replaced 
by suitable time-averages. In the case of an incompressible fluid, the thermal aspects of 
the flow are completely decoupled from the mechanical aspects, so it is legitimate to 
identify the free energy with the kinetic energy. The same is true for a compressible but 
cold fluid (i.e., T = 0), which is the case discussed in Section IV. Thus the variational 
energy principle based on the action integral W (with either 5n = 0 or Sn i=- 0) can be 
regarded as a statement of the Second Law of Thermodynamics in the special case of fluid 
flow in which the thermal aspects of the flow have either been suppressed or completely 
decoupled from the mechanical aspects. 

The variational principles presented in this section and the previous one are both the 
same. The only difference is the respective classes of variations that are admitted in the 
two cases. It should be noted that the smaller class of variations (the density-preserving 
ones of this section) give rise to the larger class of extremal flows. (The class of flows 
satisfying (64) is larger than the rectilinear flows required by (49).) The reason for this 
is that each degree of freedom that is added to the class of variations gives rise to a new 
Euler- Lagrange equation which represents a new constraint on the extremal flow. 
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VI. INTERNAL STABILITY CRITERION 

A. Convexity of Action Integral as Sufficient Condition for Internal Stability 

The remarks of this subsection will apply to the variational principle of Paper II that 
is based on the action integral £& , as well as to the principle of this paper that is based 
on W . In fact, they apply to any variational principle whose action integral can be identi- 
fied with the appropriate total time-average thermodynamic potential function for the sys- 
tem under study. In such a case, local stability of an extremal flow is guaranteed if any 
conceivable variation about the extremal trial functions causes £& to increase. (The term 
“local” refers to parameter space.) Thus, since > 0, it can be asserted that 

> 0 (67) 

is sufficient to guarantee local stability. Such a case is illustrated in Figure 3A in which 
it is assumed that Stf is parameterized in terms of the two variational parameters a and b. 
Local stability does not guarantee absolute stability. If the flow in question corresponds 
to a point in parameter space that is at the bottom of a well that is connected to a deeper 
well by a very low “mountain pass,” then sufficiently strong “jostling” of the system because 
of interaction with its surrounding momentum reservoir could cause the system to jump 
over into the deeper well. A complete statistical description would involve an ensemble 
in which each of the two extremal flows would be represented with appropriate statistical 
weights. The thermodynamical description given here and in Paper II is limited to a 
discussion of local stability of a single extremal flow. 

The other three cases shown in Figure 3 do not satisfy (67). The failure to satisfy 
this condition does not, however, guarantee that the extremal flow in question is unstable. 
The reason for this is that the flow is unstable only if a dynamically acceptable alternative 
flow with lower is available to it. This will be the case only if the “spillway” running 
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from the extremum in question leads to a lower extremum. Such a case is illustrated in 
Figure 3C. If no such neighboring extremum is available, however, a saddle-point in 
parameter space is stable “by default.” Such a case is illustrated in Figure 3B. An example 
of such a case is provided by (51) when is regarded as the total time-average free energy 
of the system. Because of the second term in the integrand of (51), the extremum can be 
a saddle point. (This is discussed in greater detail in Subsection B below.) But for the 
variational principle based on <W' the extremal trajectories are straight lines in space-time 
which are uniquely determined by the imposed terminal and boundary conditions. Thus, 
the extremum characterized by these rectilinear trajectories is unique, and so necessarily 
stable. 

Figure 3D illustrates a case in which failure to satisfy (67) (at the point P 0 ) can indicate 
neutral stability (along the line P Q - P') rather than instability. If a continuum of neighboring 
extrema (the line P G - P') exists, then because each must be tangent to a horizontal plane in 
parameter space, it follows that Stf is the same for all these extrema. Neutral wave motion is 
an illustration of such a situation. So long as the wave motion is characterized by strict equi- 
partition between its kinetic energy and its convertible (i.e., free) potential energy the total 
time-average free energy (i.e., JSO is insensitive to the amplitude of the wave motion. 

The conclusion to be drawn from these examples is that the failure to satisfy the 
condition (67) does not suffice to guarantee instability. Thus (67) is a sufficient, but not 
necessary, condition for stability. However, although instability is not guaranteed when 
fluctuations exist for which < 0, it is nevertheless true that, in the case of an 

action integral that (unlike ) includes potential as well as kinetic energy, the probability 
is very large indeed that when 5 (2) ja^< 0 there exists at least one neighboring dynamically 
acceptable extremum that has a lower value of Jaf, or one having the same value of but 
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larger kinetic energy. As to the second possibility, it should be noted that in a real fluid 
with viscosity a transition to a flow having larger kinetic energy would inevitably lead to 
one having a lower value of £# . Thus both types of neighboring flows can be regarded as 
being more stable than the original extremal flow. For this reason, it is justifiable to say 
that, when <0, the extremal flow in question is “probably unstable,” and this 

expressions will be used in which follows to indicate the violation of condition (67). 

In asserting that (67) suffices to guarantee internal stability, an extremely important 
caveat must be observed, namely the fact that this is true only if the class of admissible 
flows is limited to those satisfying the imposed terminal and boundary conditions, i.e., those 
flows for which the time-average velocity of every particle is equal to a specified value. 

This makes physical sense if the average velocity is the same for all particles of the system, 
which corresponds to a system embedded in a constant-velocity momentum reservoir. If 
the reservoir velocity has shear, however, then turbulent mixing will tend, at least over 
short distances, to eliminate this shear. The inclusion of such effects would necessitate 
enlarging the comparison class of flows involved in the stability problem. That is, instead 
of comparing t .Q^-values for the members of the single ensemble of extremal flows that 
corresponds to a given specification of the average velocity (V> t for each particle (i.e., a 
given specification of Xj and X F ), it would be necessary to consider an extended ensemble 
consisting of the sum of the various individual ensembles corresponding to different speci- 
fications of XjOs) and X F (s»). This extended ensemble will be defined so as to maintain 
the constancy of some total time-average quantity such as total flux or total angular mo- 
mentum. A simple example of such an extended ensemble will be given in Subsection C 
below. For stability with respect to the extended ensemble, it is sufficient that (67) be 
satisfied for the extended class of fluctuations that results from dropping the requirement 
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that 5X = 0 on the hypersurface bounding the integration volume (V,At), and replacing it 
by the expression for 8X0, t) on the hypersurface that defines the extended ensemble. 

In Subsection B below, the laminar and internal stability criteria will be compared 
for the case of rectilinear flow with shear, and arbitrary axisymmetric azimuthal (zonal) 
flow. In making such a comparison, it is best to regard the two criteria as corresponding 
to different modes of instability. The “internal mode” maintains the time-average velocity 
profile of the base flow (i.e., the original extremal flow), whereas the “laminar mode” leads 
to the establishment of a new velocity profile. Thus a “laminar instability” is “macro- 
scopically observable” in the sense that it either leads to growing sinuous oscillations that 
tend to eliminate the velocity shear that produced them, or else it leads to an entirely 
different stable mode, such as a convection roll, whose velocity profile is different from 
that of the base flow. An “internal instability,” on the other hand, is necessarily a fine- 
scale phenomenon because it is constrained to maintain the given average velocity profile. 
The onset of such an instability would not be “macros copically observable” in the sense 
of producing large-scale changes in the flow in the fashion of a laminar mode. Rather 
it would manifest itself in a sudden increase in the degree of turbulence present in the 
main flow. This turbulence would not be random, but would have the order of the first 
(i.e., most unstable) internal mode to be excited. If, however, the base flow were so un- 
stable that many internal modes were excited, then the turbulence would appear to be 
effectively random. 

These remarks must be tempered by the caveat that was emphasized above in intro- 
ducing the extended ensemble. The introduction of an appropriate extended ensemble is a 
way of “introducing a bit of laminar instability” into a study of the internal stability 
criterion. 
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B. 


Comparison of Internal and Laminar Stability Criteria 


1. Parallel Shearing Flow 

It follows from (48) and (67) that if only variations in kinetic energy are involved 
(inertial instability) and the extremal flow is rectilinear (X = 0), then a sufficient condition 
for internal stability is 

- (5X5X:WV) • V > 0 (68) 

throughout the space-time volume (V, At). This is a conservative criterion in the sense that 
the stabilizing effect of the term (D t 5X) z that appears in (51) has been ignored. The 
reason for this, as noted in the discussion following (51), is that this term can be made 
arbitrarily small by making At arbitrarily large. 

As noted following (67), if W is used as the action integral (thereby excluding 
potential energy), the extremal flow A = 0 is always stable regardless of the sign of the 
left side of (68). The reason for this is that, for the given terminal end boundary con- 
ditions, no alternative extremum exists. 

However, if we anticipate the results of Paper II, then it can be asserted that (68) is 
the internal stability criterion that is relevant to inertial stability of rectilinear flow of a 
compressible fluid (neglecting the possible stabilizing effects of the gravitational field). For 
parallel flow in the x -direction with velocity shear in the z-direction, (68) reduces to the 
following: The flow is stable if 

~d 2 V 

V ~ ~ < 0 (69) 

dz 2 

everywhere. Figure 4 shows examples of stable velocity profiles that satisfy the stability 
condition given in (69). Figure 5 shows examples of “probably unstable” profiles that do 
not satisfy the sufficient condition given in (69). A comparison of Figure 4A and Figure 5A 
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(A) (B) (C) 


Figure 5. Sample Velocity Profiles that Violate the Internal Stability Criterion 




shows that, whereas a jet-type profile is stable according to the internal stability criterion, 
a wake-type profile is “probable unstable.” This suggests that turbulence along the center 
line of a wake in an inviscid fluid should be significantly greater than the turbulence at 
the center of a jet whose velocity profile has equal but opposite curvature. It should be 
noted, however, that this assertion is not invariant under velocity transformations parallel 
to the flow. That is, by adding a sufficiently large counter-velocity, the wake becomes 
a jet and the jet becomes a wake. This shows that, even though perfect slip is allowed at 

r>j 

the wall, it is nevertheless necessary to regard the V in (69) as referring to the rest-frame 
of the walls. The physical reason for this will be discussed in Section VII. 

The laminar criterion that corresponds to (69) is that the absence of an inflection 
point in the velocity profile is a sufficient condition for stability. This was shown by 
Rayleigh* ^ in 1880. Tollmien*^ showed that the appearance of an inflection point in a 
symmetric profile (such as that of Figure 4A or Figure 5 A) or a boundary layer profile 
(such as Figure 4B or Figure 5C) is sufficient to guarantee laminar instability. That is, 
for the profiles involved, Rayleigh’s condition is not only sufficient, but also necessary, for 
stability. Taylor* ^ pointed out that the derivation of Rayleigh’s criterion involves imposing 
boundary conditions that amount to forbidding momentum exchange between the flow 
under study and the containing walls (or surrounding fluid in a free-boundary problem). 

The question of devising more realistic boundary conditions is the laminar-flow analog of 
the internal-stability problem of defining a suitable extended ensemble. 

The literature relevant to Rayleigh’s criterion is vast, but good reviews are available. *5-19 
(The review article of Drazin and Howard,*^ which is limited to inviscid flow, is most 
relevant to the present discussion.) 
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According to the Rayleigh-Tollmien criterion, all of the profiles shown in Figures 4 and 5 
with the exception of the one in Figure 5C are definitely stable against laminar instability 
because they have no inflection point. In comparison, the profiles of Figures 5A and 5B are 
“probably unstable” with respect to the onset of internal turbulence as indicated by the 
internal stability criterion given in (69). The two cases become unstable against laminar 
instability as well if the profiles are made more realistic by imposing no-slip conditions 
at the walls, since this introduces inflection points in the profiles. If this is done, the 
internal and laminar stability criteria coincide, which means that the onset of internal 
turbulence would coincide with the onset of the sinuous oscillations that characterize 
laminar instability in parallel shearing flow. 


For the sake of making a comparison with the case of rotational flow treated in the 
next subsection, the internal stability criterion given in (68) will now be adapted to the 
case of doubly-shearing flow shown in Figure 6, i.e., flow in the x-direction that has 
different shear gradients in the y and z directions. The left side of (68) in such a case can 
be written as follows: 


- HnM(5X6X:VVV) • V = wa(6Z— 5Y) 2 

a 

r\-/ - 

+ - (ac-b 2 )(5Y) 2 
a 


(70a) 


where 

i a 2 v i a 2 v i a 2 v 

a r-u ' ~ , b — , C — “ 

V 3z 2 V dydz V 3>- 2 


(70b, c,d) 


and w is the kinetic energy density of the extremal flow. It follows from (68) and (70) 
that the sufficient condition for internal stability of double-shearing parallel flow is 

a>0 and <.ac-b 2 )>0, (71a,b) 

For b = c = 0 this reduces to the criterion given in (69). 
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2. Axisymmetric Zonal Flow 

In this case V = V(r, z)e lfl and A = - (V 2 /r)e r where (r, i p, z) are cylindrical coordinates 
and (ej, e^,e z ) are the corresponding unit vectors. Since A ^ 0, (51) and (68) do not 
apply. In anticipation of the results of Paper II the sufficient coi dition for internal stability 
will be taken to be 

6(2) w = J / 2 nM[(6X6X:VA) - (5X5X:VVV) • V] >0 (72) 

everywhere in (V,At). (The symbol 5^ 2) w has been introduced for convenience.) As in the 

r\J ~ 

case of parallel shearing flow, the stabilizing effect of !/ 2 M(D t SX) has been neglected. The 
integrand of given in (48e) has also been omitted. The reason for this is that 

it will be shown in Paper II that this term exactly cancels a term included in the expression 
for the second variation of the potential energy. If, moreover, ri is taken to be the same 
throughout the fluid, then two terms in the expression for the second variation of the 
potential energy that depend on Vn also vanish, and the criterion for inertial internal 
stability of rotational flow of an unstratified fluid is given by (72). 


This criterion can be written in cylindrical coordinates for the case of axisymmetric 
unstratified zonal (azimuthal) flow as follows: 



1 /5R\2 

+ - (AC-B 2 ) f— J + D(5$)2 > O 


(73a) 


where 


r 2 

_ r^J 

a 2 K 

r 2 

9 2 K 

A 

9 z 2 

; B — — 

j 

K 

K 

9r9z 

r 2 

9 2 K 


r^J 

r 9K 

= 1 - — 

— • 

D = 1 - 

. — 

K 

9r 2 ’ 

K 9r 


(73b,c) 
(73d, e) 
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where K = MrV is the molar angular momentum, and (6R, 6<F, SZ) are the cylindrical com- 
ponents of the particle displacement 6X. A comparison of (70) and (73) shows the follow- 
ing correspondences: (r, <p, z ) (y, x, z ) and (A, B, C) <—* (a, b, c). The term D(5$) 2 

in (73a) has no analog in (70a). The internal stability conditions arising from (73) are 

A >0; AC - B 2 > 0; D >0 (74a, b,c) 

which include the extra condition D > 0 that has no analog in (71). 

The comparison between (70) and (73) shows that it would be conceptually preferable 
to write (70) in terms of the molar momentum MV rather than the velocity V. The linear 
and angular momentum, MV and K, are the respective constants of motion of rectilinear 
and axisymmetric flow, and the spatial gradients of these constants of motion determine 
the stability of the flow. In the rectilinear case, turbulent mixing would tend to produce 
a uniform profile (MV = constant) for which internal instability is no longer possible. 
Similarly, in the case of rotational flow, mixing would tend to produce a vortex character- 
ized by K = constant, which is just a potential vortex, i.e., one for which the velocity is 
expressible as the gradient of a scalar. (V(r) w r -1 in this case.) For K = constant, it 
follows from (73) that A = B = 0 and C = D = 1, which by (74) corresponds to internal 
stability. 

It can be shown that, for vortices for which K(r) has a simple power-law dependence 
on r and no dependence on z, the vortex is internally stable according to (74) if 

- 0.62 < (d/wK/dr) < 1.0. (75) 

For the potential and rigid vortices, dlnK/dr is equal to 0 and +2 respectively. Thus some 
subpotential vortices (for which d/rcK/dr < 0) are internally stable whereas the rigid vortex 
and all suprarigid vortices (for which d/nK/dr > 2) are internally “probably unstable.” 
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The lower limit in (75) comes from violation of the condition C > 0 whereas the 
upper limit comes from violation of D > 0. The condition on D can be made more trans- 
parent by writing it in terms of the molar kinetic energy W = V 2 MV 2 instead of the molar 
angular momentum: 


r aw 

D = — — >0. 

2W 3r 


(76) 


This shows that the vortex becomes “probably unstable” when the molar kinetic energy 
increases with increasing r. If the criterion D > 0 is suppressed by imposing the constraint 
5$ - 0 so that only the condition C> 0 remains, the upper limit in (75) becomes 1.62 
instead of 1.0. 


These results for internal stability are to be contrasted with the Rayleigh criterion for 
rotational flow ’ ^ w hich says that a sufficient condition for laminar stability of rotational 
flow against the onset of an axisymmetric perturbation (i.e., an axisymmetric convection 
roll) is that the magnitude of the molar angular momentum IKI increase with increasing r. 
Rayleigh^® derived this criterion by a simple energy argument involving the interchange of 
rings of fluid subject to the constraint that each ring retain its angular momentum. (Con- 
trast this constraint on angular momentum with a constraint on the average angular velocity, 
which is the specialization of (60) that is appropriate to the case of rotational flow.) 
Rayleigh’s argument gives a sufficient condition for stability. Synge^ * showed that this 
condition is also necessary. 

In contrasting the internal and laminar (Rayleigh) stability criteria for rotating flow, 
it should be recalled that they refer to the onset of different modes of instability internal 
turbulence as compared with axisymmetric convection rolls. Subpotential vortices are un- 
stable against the onset of axisymmetric convection rolls, but (if they are not too different 
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from a potential vortex) are stable against the onset of internal turbulence. Rigid and 
suprarigid vortices, on the other band, are stable against such convection rolls, but “probably 
unstable” according to the internal stability criterion. 

The fact that different modes of instability are involved is illustrated by the case of 
the eastward-flowing azimuthal streams of air (the westerly jet streams) in the upper trop- 
osphere. These are subject to the spontaneous onset of sinuous oscillations on a planetary 
or continental scale. Such instabilities clearly fall within the province of the laminar sta- 
bility criterion. On the other hand, the clear-air turbulence that is observed near the jet 
streams is a spontaneous onset of internal turbulence that is swept along with the main 
winds of the jet stream and has no effect on the average velocity of these winds. This 
phenomenon clearly falls within the domain of the internal stability criterion. In fact, if 
the turbulent criterion given in (73) and (74) is written in terms of spherical instead of 
cylindrical coordinates and in terms of azimuthal velocity relative to earth’s surface instead 
of the absolute molar angular momentum K, then it can be shown that the north side of 
the jet maximum (in the northern hemisphere) is the region of internal instability, whereas 
the south side is stable. Observations^ show, in fact, that clear-air turbulence occurs on 
the north side but not on the south side. This prediction of the internal criterion is to be 
contrasted with the prediction based on the Rayleigh criterion for rotational flow that the 
south side of the jet maximum is where any instability should occur^ (if it occurs at all). 

C. Flux- Averaging in Parallel Shearing Flow 

The examples of internal stability criteria discussed above refer in every case to the 
restricted ensemble of flows that corresponds to imposing the condition that 5X = 0 on 
the hypersurface bounding the space-time volume (V, At) (except for segments of rigid or 
free surfaces on which SX ♦ d 0). This restricted ensemble corresponds to the class 
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of dynamically acceptable flows for which the time-average velocity <V) t for every particle 
equals the specified value. However, the fine-scale motions that appear at the onset of 
internal instability sometimes have the effect of mixing neighboring values of <V> t and re- 
placing the individual values with a single average value. Whether or not such an averaging 
process must be taken into account is obviously a physical question that depends on the 
ratio of the characteristic amplitude of the internal motion to the distance over which a 
significant change in <V> t occurs. When the physics of the flow indicates that such mixing 
and consequent averaging of <V> t does in fact take place, then this effect can be introduced 
into the formalism by defining a suitable extended ensemble that is the sum of the individ- 
ual restricted ensembles, each corresponding to a different specification of X(sr,t) on the 
bounding hypersurface. The expression for S (1 +2) ■W' given in (48) is still valid, but now 
the terminal and surface integrals no longer vanish because 5X ¥= 0 on the hypersurface. 
Moreover, the volume integrals are also affected because if, for example, SXj = 0 and 
5X F t* 0, then choosing tj = 0 it follows that 5X # 0 for 0 < t < t F . In fact, if we specify 
5X = 5X F (t/t F ), then we simply add a constant velocity 5V = D t 5X = 5X F /t F to the 

rv-» 

original velocity V. 

The evaluation of 5 (1+2 ^W for the specification of 5X p that is illustrated in Figure 
7 (and SXj = 0) will now be carried out. This specification corresponds to the flux (or 
momentum) averaging that results from “blob-exchange” between two layers in a hori- 
zontal wind with linear vertical shear. The layers are separated by the vertical distance 
5z. The problem is normalized so that each layer contains one mole of fluid and transfers 
half of this to the other layer. Thus the total amount of fluid actually moved in the ex- 
change is one mole. The final velocities of the upper and lower layers (V+ and Vi re- 
spectively) are both equal to V Q , the original velocity half-way between them. The orig- 
inal velocities of the upper and lower layers are respectively V+ = V Q + Vi8z (dV/d z ) and 
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Figure 7. Flux-Conserving Mixing of Two Layers in a Shearing Flow. (A) Before Mixing. (B) After Mixing. 

(C) Change in Terminal Positions for Particles in Each Layer 


V_ = V Q - ViSzidV/dz). The extended ensemble in this case consists of the two extremal 

flows pictured in Figure 7 A and Figure 7B. (They are both extremal because they both 

satisfy the dynamical equation of motion, which in this case is just A = 0.) The change in 
5X on the hypersurface that changes flow A into flow B is 5Xj = 0 and 

(5X f ) + = t F e*(V;- V) = - t F e x (vaSz-j); (77a) 

(5X f )_ = t F e*(V:--V) = + t F e* (77b) 

Substitution of (77) into (48) pves the difference 6 M' Vs - W^ A . For this problem the only 
non-zero integrals are (48b) and 48 f), so 


^ = 


■/ 


MV f • 5X F dN + 


/ 


y2M(D t 6X) • (D t 5X)dNdt 


(N,At) 


(78) 


= -[(l/8)M(dV/dz) 2 (5z) 2 ]t F . 

Since the flow is steady and t F = At, it follows that 

W B -W A =(# B - # A )/At= -(l/8)M(dV/dz) 2 (5z) 2 , (79) 

a result that can also be obtained (more easily) by direct calculation of the total kinetic 
energy of flows A and B. 


This example illustrates how (48) can be used to find the difference between the 
action of a given base flow < 3¥' A that is embedded in an extended ensemble and the action 
of a particular member of the ensemble < ^¥' B that corresponds to particular values of 
5X = X B - X A on the bounding hypersurface. This difference measures the relative sta- 
bility of the two extremal flows. The minus sign in (78) means that flow B is more stable 
than flow A, which corresponds to the fact that the mixing involved in the blob exchange 
is an irreversible process. 
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If the flow of Figure 7 A occurs in a statically stable stratified atmosphere in a gravi- 
tational field, then the destabilizing effect of the wind shear has to compete with the sta- 
bilizing influence of gravity. It will be shown in Paper II that when (78), which refers to 
a transfer of a total of one mole of fluid over a vertical distance 5 z, is combined with the 
corresponding expression for the change in potential (i.e., buoyancy) energy, the resulting 
expression constitutes a sufficient condition for stability of horizontal shearing wind in a 
gravitational field that is identical to a well-known criterion that results from the laminar 
(i.e., initial-value) approach to the stability problem. ^>2 5 This illustrates the remark 
made at the end of Section VI A that use of a suitable extended ensemble provides a means 
for “introducing a bit of laminar instability” into the internal stability problem. 

Finally, it should be noted that the rectilinear flow of this example is a degenerate 
case in the sense that the averaging produced by the blob exchange may be regarded as an 
averaging of momentum, velocity, or flux. In general curvilinear flow these averages are 
all different, and it is necessary to decide which of them is appropriate to use in defining 
the extended ensemble. 



VII. REFERENCE ENERGY AND VELOCITY 


In open-ended flow problems in which V extends over an infinite range, the kinetic 
action integral W diverges. The flow in such a problem is usually of interest only in a 
limited region, as in the case of flow past an obstacle, and it would be desirable to elimin- 
ate the divergence by subtracting the known kinetic energy at large distances from the 
region of interest. This can often be done by subtracting from the kinetic energy W a 
reference energy W R defined by 

Wr = Wr(A, t) = W R [A A (x, t), t] (80) 

where the functional relation W R (A, t) is specified (but A A (x, t) is not). Thus the vari- 
ation SW R results entirely from the variation of its argument functions A A (x, t). If W R 
has the form given in (80), it can easily be shown that 

(81) 

when the boundary conditions (54) are imposed. Thus, the introduction of W R into the 
integrand of ■W has no effect on S^ 1+2)( W, i.e., 


5 U+2) 


/ 


W R nd 3 (x)dt = 0 


(V ,At) 


S< 1+2) 



W R nd 3 (x)dt 


6 Cl+2) 


L 


' * V - W r (A, t)] nd 3 (x) dt 


(v,At) 


(82) 


= 8 (1+2) W'. 


If it is possible to choose W R (A, t) so that it equals W at large distances from the region 
of interest, the integral ('W - W'r ) will have a finite value whereas W diverges, and yet 
the two integrals will have the same first and second variations. 


It is important to note that W R must be a specified function of A a and t, and not 
of x and t. A specified function of x and t would not satisfy (81). 
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Instead of introducing a reference energy W R , it is possible to introduce a reference 
velocity V R = ejV^ (A) where the Cartesian components V J R are specified functions of A a , 
but not of t. If V' is the relative velocity, then it follows from (18) that 

V' = V-V R = - (a t A A )*G A - vjjA A Gr,t)] ej . (83) 

V R is a rectilinear constant of motion, i.e., for a given particle it always has the same 
magnitude and direction (which would not be the case if the components of V R were 
defined with respect to a curvilinear coordinate system, or if = V^(A, t) instead of 
V r (A)). Because of this D t V' = D t V, so if V satisfies the equations of motion, so does 
V'. As a consequence = 5 (1) W where is the space-time integral of nW' = 

VfcnMV' •' V\ The introduction of V R does, however, change the second variation. It 
can be shown that, if the terminal and boundary conditions of (52) and (54) are imposed, 
then 

5<i+2)0f» _ 5 (i+2)^ = / y 2 M(5X6X:VVV) • V R Hd 3 (jc)dt. (84) 

* / (v, At) 

The sign of the integrand depends on the choice made for V R (A), and this is arbitrary. 

The conclusion then is that the first variation of (and hence the equation of 
motion) is invariant under a Galilean transformation of a very general kind, namely one 
that can involve a different velocity transformation for each particle. For example, V R 
could be identified with the time-average velocity <V) t , in which case W' would be the 
kinetic energy associated with the deviation of the instantaneous particle velocity from 
its time-average. The second variation, however, and hence the internal stability criterion, 
is not invariant under such a transformation, or even under the usual Galilean transformation 
for which V R is a constant vector. This means that an absolute significance can be attached 
to the second variation of only if a preferred Galilean frame is specified. This point was 
made already in Section VI B 1 in comparing the respective stabilities of the jet-type and 
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wake-type profiles shown in Figures 4A and 5 A, where it was pointed out that the preferred 
frame should be chosen to be the rest-frame of the rigid containing walls. 

This choice can be explained very simply in thermodynamic terms. Free energy is 
the energy that can be converted into useful work during the transfer of some extensive 
quantity from the system to an available reservoir of this quantity by means of appropriate 
idealized engines. In the familiar thermodynamic example, entropy is transferred to the 
coldest available entropy reservoir by means of a Carnot engine. In the cases pictured in 
Figures 4A and 5A, the extensive quantity involved is momentum, and the only momentum 
reservoir available to the fluid consists of the channel walls and everything to which they 
are rigidly attached. Thus only the fluid velocity relative to the walls is significant in 
calculating the useful work that could be extracted from the fluid by transferring its mo- 
mentum to the only available momentum reservoir, and the useful work so defined is the 
free energy whose magnitude serves as a measure of relative stability. Thus, the velocity 
of the walls plays a role analogous to that of the coldest available reservoir temperature in 
the case of heat engines. In the case of a planetary-scale analysis of the earth’s atmosphere, 
the solid massive earth is the obvious momentum reservoir, and the velocity of its center of 
mass is the appropriate reference velocity. On a local scale, however, this need not be the 
case. For example, in the case of a convection cell embedded in a large mass of constant 
velocity air (no wind shear), the air mass serves as the only available momentum reservoir, 
i.e., the only one with which the convection cell can interact. 

In summary, the internal stability criterion discussed in Section VI should be referred 
to that inertial frame in which the constant- velocity momentum reservoir available to the 
system under study appears to be a rest. In this frame the expressions as given in Section VI 
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(without the inclusion of a reference velocity) represent valid stability criteria, because 
they represent fluctuations in the free energy that could ideally be extracted from the sys- 
tem in transferring momentum to the specified momentum reservoir. 


I 


I 


I 




VnL CONCLUDING REMARKS 




A necessary preliminary to the construction of a variational energy principle for fluid 
flow is a discussion of the first and second variations of the kinetic action <9? , and of the 
terminal and boundary conditions that must be imposed on the trial functions, and this 
has been the main business of this paper. 


The total kinetic action W can serve by itself as the action integral for either an in- 
compressible fluid, or for a cold compressible fluid, both in the absence of any external 
forces. In these simple cases, 9? can be regarded as the total free-energy functional of the 
flow. In Paper II the appropriate free-energy functional that includes compressive and 
gravitational energy will be derived. 


The same terminal and boundary conditions on the trial functions describing the flow 
kinematics that have been discussed in this paper will also apply to the variational energy 
principle derived in Paper II, and the various criteria for inertial stability of the internal 
modes that were derived in this paper will also apply to the principle of Paper II. 
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FIGURE CAPTIONS 
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Space-Time Integration Volume. (A) Trajectories Cross Bounding Hyper- 
surface Only at Terminal Volumes V r and V F but not at Bounding Surface 
. (B) Trajectories Also Cross Bounding Hypersurface at Access Surface 

Dynamically Acceptable Flows. (A) Unique Solution of Differential Equations 
of Motion with Prescribed Initial Velocity Vj. (B) Ensemble of Extrema of 
Variational Principle with Prescribed Average Velocity <V> t 

Various Types of Extrema (i.e., Stationary Points P) in Parameter Space 

Sample Velocity Profiles that Satisfy the Internal Stability Criterion 

Sample Velocity Profiles that Violate the Internal Stability Criterion 

Doubly-Shearing Parallel Flow 

Flux-Conserving Mixing of Two Layers in a Shearing Flow. (A) Before 
Mixing. (B) After Mixing. (C) Change in Terminal Positions for Particles 
in Each Layer 
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